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Abstract. We give a new time-dependent definition of spline curves in R", 
which extends a recent definition of vector-valued splines introduced by Ro- 
drigues and Silva Leite for the time-independent case. Previous results are 
based on a variational approach, with lengthy arguments, which do not cover 
the non-autonomous situation. We show that the previous results are a conse- 
quence of the Pontryagin maximum principle, and are easily generalized using 
the methods of optimal control. Main result asserts that vector-valued splines 
are related to the Pontryagin extremals of a non-autonomous linear-quadratic 
optimal control problem. 



1. Introduction. 

Polynomial splines have been extensively used in several applied areas of math- 
ematics such as computer graphics and approximation theory. Since the early 90's, 
they have been used in control theory, associated to problems of aircraft control and 
path planning of mechanical systems. These applications originated the extension 
of classical spline functions to other contexts such as Riemannian manifolds. Lie 
groups, etc. 

Another line of research started with the definition of spline functions which are 
not polynomial splines. One of the first generalizations in this direction are the so 
called scalar generalized splines, which were introduced in the 50's by Ahlberg, Nil- 
son and Walsh 1 . The connection between scalar generalized splines and optimal 
control was established between 1995 and 1999. It turns out that splines are much 
more than a tool to be used in control theory. They are intrinsic to optimal control 
problems and appear naturally as minimizers of certain problems ^1 112|. 

Recently, this connection between minimality and splines was extended to a new 
class of spline functions in arbitrary dimensional Euclidean spaces jll) . This was 
accomplished by variational arguments and a more general time-invariant optimal 
control problem. Here, using tools from optimal control, we go a step further. We 
consider a class of classical linear-quadratic optimal control problems, which are 
not necessarily time- invariant, and recover, as corollaries, the previous results. 
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2. Background. 

In this section we give an account of scalar generalized splines, its connection to 
optimal control, and collect all the necessary results to be used in the sequel. 

2.1. Scalar generalized splines. Generalized splines were first introduced in the 
late 50's by Ahlberg, Nilson and Walsh Consider the linear differential operator 
of order p E N 

L = DP ■ +ap^i{t)DP-^ ■ + ■■■ + ai{t)D ■ +ao(i)- 

where each ak{t), k = 0, 1, ... ,p — 1, is a real C^-smooth function in [a, b]. The 
operator L is acting on the space C"^[a,b] of real functions defined in [a,b]. Its 
adjoint is defined by 

L* = {-IfOP ■ + {-If-^D^-^ap^t)-) + • • • - D{a,{t)-) + a^{t) ■ . 

L* is also acting on C™[a, 6] and the scalar product for which it is computed is 
given by 

{xi,X2) = / Xi{t)x2{t) dt . 
J a 

Let A : a = to < ti < . . . < t„i — b, m E N, he a partition of [a, b], fl he the family 
of real C^^'^^-smooth functions in [a,b] which are C ^''-smooth in each interval 
[U, U+i], i = 0, 1, . . . ,m - 1 and f e fl. 

Definition 2.1. The function s : [a, 6] — > R is an interpolating generalized spline 
of / associated to A and L, if s e il, s is a solution of the differential equation 
L* Lx = in each interval [t^, ti+i], i = 0, 1, . . . , m — 1, and s{t) = f{t) on A. 

Definition 2.2. An interpolating generalized spline of / is of type I if it is such 
that sW(to) = /('^H^o) and s^^\t^) = f^'^^t^), for fc = 1, 2, . . . ,p - 1. 

The function / is usually omitted from the previous definitions. Instead, one 
has to demand that, in Definition 12. II function s fulfills the interpolation condition 
s{ti) — Si, where Si, ? = 0, 1, . . . , to, are given real numbers and, in Definition 12.21 
that s fulfills the boundary conditions s^''\to) — rj^ and s'^''\tm) = Vm where 770, 
77^, k = 1, 2, ... ,p — 1, are prescribed real numbers. Then, we just say that s is 
a generalized spline of type /. The next statement collects several results about 
generalized splines of type / which can be found in ^ . 

Theorem 2.1 ( 1 ). There exists, for each set of boundary and interpolation condi- 
tions, a unique generalized spline of type I associated with the differential operator 
L and the partition A. Moreover, this generalized spline is the unique minimizer of 
the functional 

/ {Lgfdt 

J a 

among all the functions g E Q that fulfill the same boundary and interpolation 
conditions. 

Remark 2.1. There are other types of boundary conditions, described in the litera- 
ture, that also ensure the existence and uniqueness of the corresponding generalized 
spline. 
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We now give two examples for constant coefficient operators: an example of a 
cubic spline, and an example of a trigonometric spline. Let A : < 1/4 < 1 be 
the partition of the time interval [0,1]; s{to) = 3, s(ti) = 1 and s{t2) = be the 
interpolation conditions; and s(to) = ^1, ■s(^2) = 1 be the boundary conditions. We 
first consider the operator L ~ D^. The resulting spline of type I is a C ^-smooth 
function in [0, 1] such that s{t) — cu + C2it + cs^t^ + cat^ in each where 
ciz, C2i, C3i, C4i are real constants to be found. This is the classical cubic spline. 
Considering L — + 144, the resulting spline of type I is also C ^-smooth in [0, 1] 
so that s{t) — {cii + C2it) cos (12t) + {c^i + c^it) sin {12t) in each [ti, t^+i]. 




Cubic spline Trigonometric spline 

The most immediate generalization of scalar splines to curves in M" is achieved 
by simply considering vector functions g : [a, 6] — > R", the same operator L as 
before, and adapted interpolation conditions, boundary conditions, and set Q. It 
is obvious that each component of the resulting spline will be a scalar generalized 
spline, and therefore such a spline curve will always minimize the functional 




(Lg, Lg) dt , 



where (• , •) stands for the Euclidean inner product, among all functions in Q that 
fulfill the same boundary and interpolation conditions. As we shall see, from an 
optimal control perspective such a trivial generalization is not the natural way of 
extending scalar-splines to vector- valued splines. 

2.2. Scalar generalized splines and optimal control. Since the early nineties, 
in order to deal with applied problems from Robotics, there has been an increasing 
interest to combine spline curves and integral cost problems associated with linear 
control systems. Among theoretical developments, it was found that scalar gen- 
eralized splines are minimizers of a simple optimal control problem with a linear 
time-invariant control system and a single control (see [HI and |12j). This discovery 
is of crucial importance, because it introduces a new perspective to the subject: 
scalar spline functions are better viewed as a consequence of the search for an op- 
timal control, rather than a postulate imposed a priori in order to solve particular 
classes of problems. Given its importance, we summarize the main result here. 
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Consider the following autonomous linear-quadratic optimal control problem: 

f'' 

min / dt 

J a 

subject to 

x = Ax + Bu 

x(to) = Xo, x{tm) = Xm 

x^{ti) = ai, « = 1, 2, . . . ,m - 1 , 

where a = < ti < ■ ■ ■ < tm-i < t„i — b, x^ is the first component of the state 
vector, Qii e M, u is a scalar function which is C ""^-smooth in [a, 6] and C "-smooth 
in each interval Let us assume that the state space is R" and that the 

state vector is a C^"~^-smooth function in [a, b] which is also C ^"-smooth in each 
interval [ti,ti+i]. 

Theorem 2.2. // the control system x = Ax + Bu of problem is completely 
state controllable with matrices A and B in the canonical form 
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for given real numbers aj, j = 0, 1, . . . , n — 1, then the optimal control problem 
has always a unique solution with the first component of the optimal state vector 
being a generalized spline of type I associated to the constant coefficient differential 
operator L = D" — a„„iZ3"~^ _ . . . _ aiD — oq. 

If the first component of the optimal state vector of problem J^l is a scalar 
generalized spline, the following questions come immediately to our mind: What 
can be said about the minimizing state trajectory of the optimal control problem? 
Is it some sort of a generalized spline in R"? The answer to these questions leads 
us (see Definition 13. 1|) to a new time-dependent definition of generalized spline in 
M", which is the main contribution of the present paper. 

2.3. Pontryagin's niELximum principle, existence, and regularity. The gen- 
eral problem of optimal control can be defined, in Lagrange form, as follows: 

, min I[x{-),u{-)]= r C{t,x{t),u{t)) dt 

(2) x{t)^^{t,x{t),u{t)) 

{x{ta),x{tb)) = (q;,/3) 

a;(-)eM^i,i([ia,ifc];R") 
GZ^([t„ifc]; riCM'-) . 

We assume that C : [ta,tb] x W x W ^ R and ip : x R" x R*^ ^ R" 

are C^-smooth functions with respect to all arguments, and that the boundary 
conditions, together with the class of control functions U, are given. The standard 
method to solve Q is usually based on the deductive approach: (i) a solution 
exists for the problem; (ii) the necessary conditions are applicable, and they identify 
certain candidates (so called extremals); (iii) subsequent elimination (if necessary) 
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identifies the solution (or solutions). We are interested in the case where there 
are no restrictions on the control variables: = M'". The unrestricted case poses 
many difRculties, and the problem turns out to be a difficult one, even in special 
situations. As we explain next, most part of difficulties appear in the application 
of steps (i) and (ii). 

The first general answer to (i) was given by A. F. Filippov in 1959 7 , assuming 
the admissible controls to be integrable {U — Li), and the control set H. to be 
compact. As far as we assume to be a noncompact set, Filippov's theorem does 
not apply. To solve the existence problem, we make use of the following theorem 
(see m). 

Theorem 2.3 ("Tonelli" Existence Theorem for Q). Problem Q has a minimizer 
{x{-),u{-)) with u{-) € Li {[ta,th]: W), provided there exists at least one admissible 
pair, and the following convexity and coercivity conditions hold: 

• (convexity) Functions C(t, x, •) and ip{t, x, •) are convex for all {t, x); 

• (coercivity) There exists a function 9 : M, bounded below, such that 

C{t,x,u) > 6 {\\(p{t,x,u)\\) for all {t,x,u); 
fim = +oo; 

r-^+oo r 

lim \\ip{t,x,u)\\ = +00 forall{t,x). 

|t«lH+oo 

Remark 2.2. For the definition of convexity of C{t, x, •) and (p{t, x, •) see 0]. In the 
case if = u one has the fundamental problem of the calculus of variations, and we 
get from Theorem 12 .31 the classical Tonelli existence theorem. 

Step (ii) is addressed by the Pontryagin Maximum Principle 

Theorem 2.4 (Pontryagin Maximum Principle). If {x{-), u{-)) is a minimizer of 
^ and u(-) is essentially bounded, u(-) € L^o, then there exists (ipg, tpi')) 7^ 0, 
"00 5: 0, V'(') S such that the quadruple (x(-), u(-), V'Oj '/'(■)) ^■^ Pontryagin 

extremal: it satisfies 

• the Hamiltonian system 

(3) 



_ dH 

^ _ dH 
dx 



the maximality condition 



(4) H [t,x{t),u{t),^o.i^{t)) ^ ui&^H {t,x{t),v,^o.i^{t)) ■ 
with the Hamiltonian 

(5) H{t, X, u, ipo, ip) ^ipoC {t, X, u) + {tp, if (t, X, u)) . 

Definition 2.3. A Pontryagin extremal {x{-), u(-), i{jq,4'{')) is said to be abnormal 
when ipo is equal to zero, and normal otherwise. 

The existence is assured in the class of integrable controls {U — Li), while the 
formulation of the Pontryagin maximum principle assume the optimal controls to 
be essentially bounded (U ~ Loo C Li). For minimizers predicted by existence 
theory. Theorem 12.41 may fail to be valid, because the values of optimal controls 
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can be unbounded. This is a possibility even for very simple instances of problem 
e.g. L a polynomial and 93 linear. One such example can be found in the 
problem 



dt 



(6) 



x{t)^u (t) 
a;(0) = 0, x{l)^k, 



satisfies all the hypotheses of Theorem l2.3l it can be proved (see that for certain 
choices of constants k and e there exists a unique optimal control u (t) = kt^^/^\ 
but Theorem 12.41 (Pontryagin maximum principle) is not satisfied since Tp{t) — 
Cx {t,x{t) ,x{t)) = ct^"^/^ is not integrable (V'(-) is not an absolutely continuous 
function). 

In order to apply the deductive method (i)-(iii) one needs to close the gap 
between the hypotheses of existence and necessary optimality conditions. For that, 
conditions beyond those of convexity and coercivity, assuring solutions u{-) to be 
in Lao and not only in Li, must apply. To exclude the possibility of bad behavior 
that occurs for ©, we will focus our attention to problem ^ with Lagrangian £ 
and function ip given by 

£{t,x,u) = {B{t)u,B{t)u), 
^ ' ip{t,x,u) = A{t)x + B{t)u, 

under the hypothesis 

(HI): B{t) is a square matrix with full rank for all t; 

(H2): the dynamical control system x{t) = A{t)x{t) + B{t)u{t) is completely 

state controllable; 
(H3): t — > A{t) and t B{t) are C^-smooth functions. 

Roughly speaking, this gives the biggest class of optimal control problems which 
generalize in a natural way; do not admit abnormal extremals (see the next re- 
mark) ; and for which the gap between existence and necessary optimality conditions 
is automatically closed. 

Remark 2.3. Since there is no constraint on the control, singular trajectories are 
exactly projections of abnormal extremals. But due to the assumption on the linear 
system (it is supposed to be completely state controllable), there is no singular 
trajectory, and thus the optimal control problem has no abnormal extremals. 

Theorem 2.5 (Boundedness of optimal controls ^Sl)- Under the hypotheses of the 
existence Theorem \2.S[ if there exist constants c > and k such that 



dL 



dt 



dt 



< c\C\ + k. 



<cM+k, 



dC 



dx 

dip, 



dx 



< c\C\ + k, 
<c\ipi\+k (i = 1, 



n). 



then all 



mimmizers 



of ^ satisfy the Pontryagin maximum principle. 



It is a simple exercise to see that with C and defined by Q, hypotheses (HI) 
and (H3) imply all the conditions of Theorems 12 . 31 and 12 . 51 
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3. Main results. 

We are interested in the following non-autonomous linear-quadratic optimal con- 
trol problem: 

min J[u{-)]^ [ {B{t)u{t),B{t)u{t))dt 

(p) subject to 

x{t) = A{t)x{t) + B{t)u{t) 
x{ti) = x^, z = 0, l,...,m, 

for a given partition a ^ to < ti < ■ ■ ■ < tm-i < tm — b and fixed Xi £ M". 
The control u : [a, b] — > M" is unrestricted; the state function x : [a, b] M" is 
an absolutely continuous function; A{t) and B{t) are n x n matrices and B{t) is 
nonsingular. We find the minimizer of (P) by solving (Pi), i = 0, 1, . . . , to — 1, in 
the interval 

/ti+i 
{B{t)u{t),B{t)u{t)) dt 

(P) 

^ x{t) = A{t)x{t) + B{t)u{t) 



x{ti)^Xi, x{ti+i)=Xi+ 



In order to guarantee the applicability of the Pontryagin maximum principle, and 
the existence of a normal solution, hypotheses (HI), (H2) and (H3) of previous 
section are in force. Under these assumptions we can choose, without any loss of 
generality, ^/^o = — 5 in Theorem 12.41 The Hamiltonian (jSJ is then given by 

H{t, x, u, ip) = -\u'B{t)'B{t)u + iy{A{t)x + B{t)u) , 

where we use the symbol prime ' to denote the transpose of a given vector or matrix. 
The Hamiltonian system Q reduces to 



(8) 



j x{t) = A{t)x{t) + B{t)u{t) , 

\m^~A{tym, 

while from the maximality condition Q one obtains 

B{t)'{B{t)u{t) - i;{t)) =0. 
This equation implies that V-'(i) — B{t)u(t) and hence 

(9) u{t) = B{t)-^m 

is the unique Pontryagin extremal control. Thus, due to Theorem 12.51 u given by 
(jSJ must be optimal. 

From equation ■4>{t) ~ B{t)u{t) and from equation \lj{t) = —A{tyip{t) of system 
^ we get the matrix differential equation 

(10) j^iBit)uit))+A{tyBit)uit) = 0. 

Introducing the matrix differential operator L = D — A(t), the control system 
x{t) — A{t)x{t) + B{t)u{t) can be written as 

(11) Lx{t) = B{t)u{t) 
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and equation (fTn|l as 

(12) L*B{t)u{t) =0, 

where L* — —D — A{ty is the adjoint operator of L. From (|11() and (|12|l we conclude 
that the minimizing state trajectory is a solution of the differential equation 

L*Lx{t) = 0, 

which can be written as 

£{t) + {Aity - A{t)) x{t) - {A{tyA{t) + A{t)) x{t) = . 
We have just proved Lemma l3. II 

Lemma 3.1. Under hypotheses (H1)-(H3) the optimal control u is, in each interval 
[ti, ti^i], i = 0, 1, . . . , m—1, a solution of the matrix differential equation L*B{t)u = 
with L* = —D~A{ty the adjoint operator associated to the operator L = D — A{t). 
The corresponding optimal state trajectory x is such that L* Lx — in each interval 

[ti,ti+l] ■ 

An explicit expression for the optimal state trajectory and for the optimal control 
can be obtained in terms of the state transition matrix. These results are stated in 
the following Theorem. We refer the reader to ,3., for the definition, and properties, 
of the state transition matrix. 

Theorem 3.2. The optimal state trajectory of problem (P) has, in each interval 
[ti, ti+i], i ~ 0,1, . . . ,m — I, the explicit expression 

(13) x{t) = ^it,t,)x, + (^j^ $(t,s)$(t,,s)'ds^ S-\<^{t„t,+i)x,+i - x,) 

where $ is the state transition matrix associated to x = A{t)x, and S is the sym- 
metric matrix given by 

<^{ti, s)<^{ti, sy ds . 



Furthermore, the optimal control of problem (P) has, in each interval [ti,ti^i], 
i = 0, 1, . . . , m — 1, the explicit expression 

(14) u{t) = B{t)-^ $(i„ t)' S-\<i>it,,t,+i)x^+i - Xi) . 

Proof. fTheorem l3.2|l Since -0 = B{t)u, the Hamiltonian system takes the form 



(15) 



X = A(t)x + -0 , 

tP = -A{ty^. 



From equation = —A{tyilj we get tp{t) = $(ii, ty'ip{ti). The substitution of ip in 
equation x — A{t)x + if) of system (jlSII generates x = A{t)x + ^{ti,tyil){ti). The 
solution of this complete differential equation, with initial condition x{ti) — Xi, is 
given by 

(16) x{t) = $(t, t{)x, + ^{t, s)$(t„ s)V(i») ds. 
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Now, we just have to find ipiti). Using the other initial condition x{ti+i) — Xi+i 
we get 

Xi+i = ^{U+i,ti)xi + (^J^ $(ii+i,s)$(ii,s)'(is^ tPIU) 



ti 



If we denote the symmetric matrix 

' ^{u,s)^{t,,sy ds 



by S'(ti,ti+i), or simply by S, we can write 

$(ti+i,tj)~^a;j+i -Xi = S ij){ti) ^ <I>(ij,ti+i)a;i+i - Xi = S ipiU). 

Since matrix S is always non-singular, wc get ip{ti) = S^^{^{ti,ti^i)xi-i-i — Xi). 
Finally, from H16(l . we obtain the equality Ijl^ll : 

a^(i)|tG[t,,t,+i] = U)x^ + s)^{U, s)' ds^ tj+i)x,+i - a;,). 

The second part of the theorem is a direct consequence of equation J^. From 
previous calculations we have 

ij{t) ^ t)' S-\<^(U, U+i)x,+i - X,) 

and thus, equality (|14(l follows immediately. □ 

Remark 3.1. From the proof of Theorem 13. 21 it follows, by direct calculations, that 
the optimal value for the integral functional Ji[-] of problem {Pi) is given by 

{^{ti,ti+i)xi+i - Xi)' S^^ {^{ti,ti+i)xi+i - Xi) . 

Remark 3.2. We have seen that in each interval [ti, i^+i], i = 0, 1, . . . , to — 1, the 
optimal state trajectory of problem [P) is a solution of the matrix differential 
equation 

x{t) + {A{ty - A{t)) ±{t) - {A{tyA{t) + A{t)) x{t) = 

which does not depend on the matrix B{t). This is natural since we can make the 
substitution v = B(t)u in the problem (P) and thus eliminate the presence of 
matrix B(t) in all further calculations. 

Remark 3.3. When problem (P) is autonomous, the first part of Theorem 13.21 
reduces to Theorem 2.12 in 

Lemma 13.11 and Theorem 13.21 give the main motivation for our definition of 
generalized time-dependent spline in R". Let L be the linear matrix differential 
operator of order p 

(17) L = DP-- Ap^,{t)DP-^ A,{t)D ■ -Ao{ty , 

where each Aj{t), j = 0, 1, . . . ,p — 1, is a real square n x n C^-smooth matrix 
function in [a,b]. The operator L is acting on the space C""[a,6] of real vector 
functions defined in [a,b]. The adjoint of L, denoted by L* , is defined as 

L* = {-1)PDP ■ + {-l)PDP-\A;^,it)-) + (-l)P-^DP-\A;_,{ty) + 

+ ■ ■ ■ + DiA[ity) - A',it) ■ . 
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L* is also acting on C "'[a, 6] and the scalar product for which it is computed is 
given by 

{xi,X2) — I Xi(t)' X2{t) dt . 
J a 

Consider 

(18) A : a = tQ < ti < . . . <tm ^ b 

to be a partition of [a,b], and let represent the set of all R" -valued functions 
defined in [a,b] which are of class C^p~^ in [a,b] and of class C^^ in each interval 

[U,U+l], Z = 0, 1, . . . ,TO - 1. 

Definition 3.1 (Generalized time-dependent spline in K"). A function s : [a,b] 
K" is an interpolating generalized spline of / G fi, associated to A (I18|) and L 
(|17|l . if s € ri, s is a solution of the matrix differential equation L*Lx ~ in 
each interval s{t) — f{t) on A (interpolation conditions), and s''^''(io) = 

/W(*o), s^'^^tm) = f'^^\tm), for A: = 1,2, ... 1 (boundary conditions). 

Remark 3.4. Definition 13.11 includes, as particular cases, the scalar Definition 12.21 
and the definition introduced in lli . 

Remark 3.5. As done in the scalar case, the interpolating function / G 51 can be 
omitted in Definition 13. II 

Remark 3.6. The function x{t), t G [a,b], given in each interval [tj,ti+i], i — 
0,1,...,TO — 1, by H13|l . is a generalized time-dependent spline in M" in the sense 
of Definition 13.11 

Remark 3.7. For L = the solutions of L*Lx = give polynomial splines in R" 
with all the components being scalar polynomial splines of degree 2p — 1. This is, 
as mentioned at the end of t |2.1l the immediate generalization of scalar polynomial 
splines to vector-valued splines, and the one found in the literature. 

We have seen that generalized splines associated to an operator L of order p = 1 
are related to the optimal control problem (P). For p > 1, there corresponds an 
optimal control problem with higher-order dynamic x^^^ = X]^=o ^ji^)^^''^ +B{t)u. 
This higher-order optimal control problem can be easily written in form (P). For 
that we introduce new state variables, reducing the control system of order p to a 
first-order control system. This is the same to say that when L is an operator of 
order p > 1, the homogeneous differential equation L* Lx = of order 2p can be 
reduced to a first order differential equation, just by increasing the dimension of 
the matrices Aj{t), j = 0, 1, . . . ,p — 1. 

Under our hypotheses, it is possible to write the optimal control problem (P) 
as a problem of the calculus of variations with higher-order derivatives. This is 
done by showing that an arbitrary admissible pair {x{-),u(-)) of (P) can be always 
expressed in terms of higher order derivatives of a single vector valued function (see 
[S]). From Theorem 13. 21 we obtain: 

Theorem 3.3. Given the operator L (|17|) and the partition A (|18ll . there exists a 
unique generalized spline in R" for each set of boundary and interpolation condi- 
tions. This generalized spline is the unique solution of the following higher-order 
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problem of the calculus of variations: 

rb 

{Lg,Lg)dt — > min , 

among all the functions 5 G that satisfy the same boundary and interpolation 
conditions. 

4. Examples. 

We give two examples for which the state and control spaces are . We denote 
the components of the state vector x by xi and X2; the components of the control 
vector u by ui and U2. The first example is 

2 

2 



min / ui + U2 dt 

tl = (Ml,M2)' Jo 

subject to the control system 

±1 = t^X2 + U2 , 



(19) 



X2 — —t^Xi + Ml 



and the interpolating conditions 

x(to = 0) = (0,0)', a;(ti = f) = (f,0.5)', a;(t2 = 2) = (-0.25, 1)' . 

As far as the control system ()19|l is non-autonomous, this example is not covered 
by the results in We have the time interval [0, 2] and its partition A : a = < 
1 < 2 = &. The associated state transition matrix is given by 



$(t,i,) 



cos 



(^) 



sm I *' *' I cos 



The linear dynamic is completely state controllable. Such a conclusion follows 
immediately from the fact that the symmetric matrix 

r-Ti 



W= j \{T^,s)B{s)B{sy<^{TQ,sy ds 



is positive definite for some ti > tq with to,ti S [0,2]. This is a classical test for 
complete controllability which is due to Kalman |S] . Since B and $ are orthogonal 
matrices, the matrix W is simply 

n - To 

Tl - To 

The optimal control is, in each interval [ti,ti-|_i], solution of the equation 

L*Bu = ^ 

We get 

sin cu + cos (^-Y^) C2t 



ill + t'^U2 = 

U2 — t^Ui = 0. 



(20) u{t)\t^[t^^t.^^] = 



— sm I 

cos (^-g^jcii + sin (^^-^^C2i 
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where cu and C2i are real constants to be found. The corresponding generahzed 
sphnc in , solution of equation 

r*r — n<t^/^i~ 2f^i2 - t^xi - 2tx2 = 

\ X2 + 2e±i + 2tXl - t*X2 = 

in each interval [U,ti+i], is given by x{t)\te[t,,t,+i] = (xiit), X2(t)y with 

Xi{t) = cos (^-3^) (xiiU) + {t~ U) Cii) + sin (^^^) {X2{ti) + (i - U) C2i) 



and 



X2 



(t) = - sin (^^) (xi {U) + {t- U) cu) + cos (^^) (2:2 (i,) + (t - t,) C2i) 



where cu and C2i are the same constants that appear in formula H20|l . As expected, 
the resulting spline is a continuous vector function and the optimal control function 
is discontinuous at t = ti. 




First example - optimal control 



GENERALIZED SPLINES IN K" AND OPTIMAL CONTROL 
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We now apply our results to the autonomous situation treated in (IJj. Consider 
the optimal control problem 

c4 



min I ui + 2uiU2 + 2u2 dt 



subject to 

r il = -X2 + U2 

[ ±2 = 2X1 +Ui+ U2, 

x(to = 0) = (0,0)', xih = 1) = (1,0.5)', 

a;(t2 = 2) = (-0.25,1)' and a;(t3 = 4) = (1, -1)'. 

The optimal state trajectory is the generalized spline which, in each interval [U, ti+i], 
is solution of equation 

L*Lx ^O^x+{A' - A)±~ {A' A) a; = ^ ( ^Ji + - 4x1 = 

[ X2 - 3X1 - X2 = 0. 

We get, in each interval [ti, i^+i], 

Xi{t) = sin (\/2t)(f Ci, + Cu) + cos {y/2t){^Cu + f C2» + C3^) 

and 

X2{t) = sin iV2t){cu + ^C2^ + V2C3^) + COS (V2<)(-^ci, + ic2, + %/2 C4O 

where cu, C2i, c^i and ca are real constants to be found in each interval. 



Second example - generalized spline in 
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